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Abstract. 

C\ , We consider solutions to the hyperbolic system of equations of ideal granular 

.^ | hydrodynamics with conserved mass, total energy and finite momentum of inertia 

and prove that these solutions generically lose the initial smoothness within a finite 
time in any space dimension n for the adiabatic index 7 < 1 + §• Further, in the 
one-dimensional case we introduce a solution depending only on the spatial coordinate 
outside of a ball containing the origin and prove that this solution under rather general 
assumptions on initial data cannot be global in time too. Then we construct an 
exact axially symmetric solution with separable time and space variables having a 
strong singularity in the density component beginning from the initial moment of 
\n ' time, whereas other components of solution are initially continuous. 
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1. Introduction 

The motion of the dilute gas where the characteristic hydrodynamic length scale of the 
flow is sufficiently large and the viscous and heat conduction terms can be neglected is 
governed by the systems of equations of ideal granular hydrodynamics [l]. 
This system is given in R x R", n > 1, and has the following form: 

d tP + dw x {pu) = 0, (1.1) 

d t (pu) + Div x (pu ® u) = -V x p, (1.2) 

d t T + (u, W X T) + ( 7 - l)Tdiv x u = -ApT 3 / 2 , (1.3) 

where p is the gas density, u = (u\, ..., u n ) is the velocity, T is the temperature, p = pT 
is the pressure, and 7 is the adiabatic index (1 < 7 < 1 + |), A = const > 0. we 
denote Div^ and div x the divergency of tensor and vector, respectively, with respect to 
the space variables. The only difference between equations ( ll.ip -f TOj) and the standard 
ideal gas dynamic equations (where the elastic colliding of particles is supposed) is the 
presence of the inelastic energy loss term —ApT 3 / 2 in ( II. 3p . 

The granular gases are now popular subject of experimental, numerical and 
theoretical investigation (e.g. [T], [7], [I] and references therein). In contrast to 
ordinary molecular gases, granular gases cool spontaneously because of the inelastic 
collisions between the particles. The inelasticity of the collisions generally causes the 
granular gas to form dense clusters. The formation of complex structure of clusters 
has been investigated by means of molecular dynamics simulations and hydrodynamic 
simulations. 

The Navier-Stokes granular hydrodynamics is the natural language for a theoretical 
description of granular macroscopic flows. A characteristic feature of time-dependent 
solutions of the continuum equations is a formation of finite-time singularities: the 
density blowup signals the formation of close-packed clusters. 

System ( II. ip - ( 11.31) can be written in a hyperbolic symmetric form in variables 
p, u,K — pp' 1 and therefore the Cauchy problem 



p,u,pp-^\ t=0 eH m (W l ),m> 1 + 



n 
L2J 



has a solution as smooth as initial data at least for small t > [6 J . 

We will call the solution to (jl.ip - (I1.3P classical if p > 0, p > and the components 
of solution belongs to C\[0,T), H m (R n )), T < 00. 

System ( II. ip - ( II. 3p has no constant solution except the trivial one (p = 0), therefore 
the solution with the components highly decreasing as |x| — > 00 can be considered as a 
natural perturbation if this steady state in the case of the mass conservation. 

Let us note that there exists a solution (the homogeneous cooling state) with 
constant p, u and p 7^ 0. In this case the temperature T = T(t) = (-^ + T(0) -1 ) -2 , 
where T(0) is the initial value of temperature (the Haff's law). Another trivial solution 
is u = p = T = 0, pit, x) = po(x). 
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We introduce the following integrals: the total mass 

M{t) = pdx, 

the momentum 



jw = /,»*, 



and the total energy 



f(t) = J ( ^p\u\ 2 + pT) d.r = E,.(t) + E-M). 

Here E^(t) and Ei(t) are the kinetic and internal components of energy, respectively. 
Let us introduce also the functionals 

G(t) — 7: p(t, x)\x\ 2 dx, F{i) = I (u,x)pdx, 

where the first one is the momentum of inertia. 

We consider below the solutions to (II. ip - (11.31) such that the integrals m, P, E and 
G converge and call them solutions with finite momentum of inertia (FMI). It is easy 
to verify that for this class of solutions M{t) = M = const, P(t) = P = const, 

E'(t) = — [pWpWdx. (1.4) 

J-i- J 

R n 

The latter equation expresses the inelastic energy loss per collision. 

Further we consider a function K = pp~ 7 , where \nK can be interpreted as the 
usual hydrodynamic entropy. System (II .ip - (11.31) result 

f| = _AA-§p¥< , (1.5) 

therefore 

K(t, x) < K + = sup K(0, x). (1.6) 

2. Main theorem: nonexistence of global smooth solutions 

Theorem 1 Let P^O and M be sufficiently small. Then there exists no global in time 
classical FMI solution to the Cauchy problem for U.l\) - U.3\) . 

To prove the theorem we need to get firstly certain estimates of energy. 

Lemma 1 For the classical FMI solutions to U.l\) - U.3\) the following estimates hold: 

P 2 

E k (t) > -— = const, (2.1) 

v ; ~ 2M 

37-1 



Z\t)<-hC x Ef^{t), (2.2) 

where C x = K+~ (7 - l)^iT M" 5 ^!. 
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Proof. Inequality ( 12.1 j) follows immediately from the Holder inequality. To prove 
(12. 2p we firstly use the Jensen inequality as follows: 

/ pdx\ ^FTT) / j Kpp~ / ~ 1 dx\ ^^ 



M J \ M 

J pK*^~) p^r~ dx J p x l 2 p z l 2 K^ 1 dx -^ 

< ^ < *: < ^±_ f P V2 p 3 /2 dx _ (2 3) 

M ~ M ~ M J H F y ' 

R" 

Together with (O) inequality O gives (E2]). D 

The following lemmas establish the properties of the momentum of inertia. Acting 
as in [2] we get 

Lemma 2 For classical FMI solutions to U.l\) - U.3\) the equalities 

G'(t) = F(t), (2.4) 

F'(t) = 2E k (t)+n( 1 -l)E t (t). (2.5) 

take place. 

Proof. The lemma can be proved by direct calculation using the general Stokes 
formula. □ 

Then we get two-sided estimates of G(t). 

Lemma 3 Ifj < 1+-, then for the classical FMI solutions to U.l\) - U.3\) the estimates 
P 2 



2M 
hold. 



r + F(0)t + G(0) < G(t) < 8(0)1? + F(0)t + G(0) (2.6) 



Proof. First of all (EU, (EZJ) result 

G"(t) = 2E k (t) + n( 7 - l)Ei{t) = 28(1) - (2 - n( 7 - l))^(t). (2.7) 

Therefore together with (j2.ip we have 

^ < G"(t) < 2£(0), (2.8) 

after integration this gives (I2.6p . D 

Now we get an upper estimate of Ei(t). 

Lemma 4 If 7 < 1 -+- - and ?^ 0, t/ien /or the classical FMI solutions the following 
estimate is true: 

E W < g£=w> (2-9) 
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Proof. The method of obtaining the upper estimate of Ei (t) is similar to [2j • Namely, 
let us consider the function Q(t) = 4G(t)S(t) — F 2 (t). The Holder inequality gives 
F 2 < 4G(t)E k {t), therefore £(t) = E k (t) + E^t) > E^t) + gj| and 

W) < ||. (2.10, 

We notice also that Q(t) > provided the pressure does not equal to zero identically. 
Then taking into account ( 12.41) . ( 12.51) and (1 2. 7ft we have 

Q'(t) = AG'(t)S(t) - 2G\t)G"{t) + AG{t)£\t) 

= 2(2 - 71(7 - 1)) G'{t)Ei{t) + AG(t)£'(t). (2.11) 

Further, one can see from ( 12 .8p that G'{t) > beginning from a positive to for all initial 
data. Thus, for 7 < 1 + f we have from (12301) . (12111 

qm 2-n(7-l) gff) 

g(t) - 2 c(t) • ^ J 

Then d2~T0j) and ( 1212]) give 

^ W -G(7-iW2(i)' C2 " 4 ' 

The proof is over. D 

Remark 1 Inequalities Ii2. 1\) and Ii2.2\) result 

37-1 
£ '( f )<-A Cl (£ W -_^ 

Integrating this inequality we the following upper estimate of Ei(t) : 

E i (t)<( Cl t + c 2 y 2J ^ 1 , 

where c\ = — ^ — - and c 2 = (E(0) — P 2 /{2M)) T+ 1 are positive constants. However, 
for 7 < 1 + - this estimate is less exact than Ii2.9\) and is not enough for our proof. 

The next step is a lower estimate of Ei(t). 

Lemma 5 Let P^O. Then for the classical FMI solutions the estimate 

holds with a positive constant 

1 7(n + 2)-2 ( 7 -l) 7 (n + 2)-2 

the value of Cq is written in /i2.19\) . 
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Proof. The proof is based on the inequality 

27 "(7-1) 



11/11* 


(It' 

( 


2 7 




.,11/11 

n(-y-l) 


(n + 2)7 
LT(M": 

■71 

+ 1 


X) 11/11 
/ 2 7 


(n + 2)7 

L!(K"; 

V 


— n 
|x| 2 dx)' 

-27 


c 


(n + 2) 7 - 


,n + 2)7 — n 


7l(7- 


dJ 


Wt~ 


l)J 





established in [2]. 

Namely, we have for / = Kp 



EAt) = — ?— ■ [ 2- ^- dx > ^ / (KpV dx 

w 7 - 1 J p 7 K-y - k\- x (7-1) ; • '■ 



> 



7(u+2)-2 

fori r Kodx^l 



R™ / 



^r 1 (7 - 1) / \ t 



7 (n+2)-2 

> I (SJ£W) ' (214) 

- 2 h-iW2A'<?'" 1)(1+ " /2) (7-l) (G^)) 1 ^ ' 



where we denoted £?(£) = f Kpdx. 

R" 

Further, from ( II. ip and (II. 5p we have 



therefore 



d t (Kp) + div x (Kpu) = -Aid /r* 2 , 



S'(t)= -A / K'lp^dx. (2.15) 



From the Jensen inequality we get 



1+7 



3 / f Kp~<dx\ ^^ 

37+3 KIM f 7+1 3 ' X 

K2 p— dx < + / K pp— dx < K\M 



S(t) J rr " + 5(t) 



1+7 
) 



/^(t)y(^ 



(2.16) 



where C 4 = #1 M (7 - 1)^) . Further, from (|23j) . ( I2TT5J) . (1216]) we obtain 

1+7 7 + 1 nf-Y+1) 

S*&=V(t)S'{t) > -AdC^ {G{t))~*. (2.17) 

Now we take into account the lower estimate in (12. 6p together with the fact that 
beginning from a certain to the value of F(t) becomes positive if P 7^ (see (12. ip . 
( 12. 5p ) and integrate (I2.17p . Thus for t > t we get in the case 7 < 3 

2(7-1) 



7-3 4(7 - 1) 



,1- 



n(7 + l) \ 7-3 



W)*[m*™ + K ^J C,- • ] (2.18) 
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where C 5 = M (F(* ))~^~^ P~ 2 , and for 7 = 3 {n = 1) 

Sit) > S{t ) exp f^\ . (2.19) 

We denote the constant in the right-hand side of (I2.18p . (12. 19[) by C^j). 

It is easy to see that for sufficiently large t the value of S(t) is separated from zero. 
Thus, from fl2TTij) . (12181) . fl2TT9|) we obtain f l2TT3|) . □ 

Proof of Theorem^ Taking into account (12. 2p . fl 2 . 1 3 [) and the right-hand side of (12.61) 
we get 

S\t) < -Ad cfi^ (Git))- 111 ^ 1 (2.20) 

< -ACi cfi^ (S(0)t 2 + F(0)t + G(0))- :iI ^ ii . 

As follows from (12.11) . ( 12.41) beginning from to > the value of F(t) becomes positive. 
Integrating (I2.20J) we obtain for t > t 

37-1 2 n ( 3 7-l)/2('F('+ n ^-n(37-l)/2+l ( 'C ( 'Q^n(7-l)/4-l 

5(t)<5(*o)-AC7 1 C73 ,FFiy - n (37 -i)-2 + A(t) ' 

where A(t) = 0(t _n( - 37_1 ' ) / 2+1 ), t — > 00. Since C3 tends to a positive constant as M — > 
and Ci contains M in a negative degree (see the statements of Lemmas [1] and [51 then 
choosing M sufficiently small, we can always get a contradiction with inequality (12. II) . 
This prove the theorem. □ 

Remark 2 Main idea of this paper can be found in J2|/ ; where it was proved nonexistence 
of global smooth solutions to the compressible Navier-Stokes equations. For the freely 
cooling gas, it would be possible to add viscous terms (with constant viscosity coefficients) 
to apply the same technique and obtain the analogous nonexistence result. 

3. One-dimensional case 

As it was noticed system (II. ip - (I1.3P has no constant solution except of the trivial 
one p = 0. However, it is possible to construct the nontrivial steady state solution 
p(0,x),v(0,x),p(0,x) for the regions \x\ > R(t) > 0. For \x\ < -R(O) we chose the 
functions p(0, 2), v (0, x), p(0, x) arbitrarily to get initial data, smooth on the whole real 
axis. We are going to show that such solution necessarily loses its initial smoothness. 

3.1. Automodel solution 

Let us find a solution that depends on the automodel variable ^ = x — at, a = const. 
The continuity equation (II. ip gives the connection between the velocity and density as 
follows: 

w(0 = —777 + a, c\ = const 7^ 0. (3.1) 
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Equation (jl.2p and ( 13.1 P results 

ciw(0 +K0 = c 2, c 2 = const. (3.2) 

From ( 13. ip and ( 13. 2ft we have cf + p» = p(c2 — aci), therefore c 2 — ac\ > 0. 

Further, we substitute the functions u and p, found from (3.1) and (3.2) and 
expressed through p, in the equation 

d t p + ud x p + ^pd x u = —Ap 2 p 3 ' 2 , 

which is a corollary of ( 11.11) . ( 11.31) and the state equation p = pT. Thus we get an 
ordinary differential equation 

, (n A P 2 (0 ((c 2 - acMO ~ c\f' 2 

PlU " d c?( 7 + l)- 7 (c 2 -a Cl )p(0 ' 

or 

A /»<?<««> + ■*>' , (3.3) 

ci(c 2 -oci) q-7 2 (?) 

where z = (c 2 — aci)p — cf . The case c 2 — aci = that seems simpler corresponds to a 
negative pressure and we do not consider it. Equation ( 13. 3 j) has a solution 



^2 



oci / (7 + 3) arctan *I (7 + 1)V^ 2 
Ci A I ci z + cf a/z 



e = /W:= : VP " ~+ " ' T + J 7z\ + c 3, (3-4) 



C3 = const. One can see from (13. 3p that if z(£o) <£* = —, then there exists a branch 
of solution (13 .4p defined on the semi-axis (£o> +00) for positive C\ where the function 
z(£) decrease monotonically from the value z(£o) to zero. For negative ci this branch is 
defined on the semi-axis (— oo,£o), where the function z(£) increase monotonically from 
zero to z(£o)- 

3.2. "Finally steady state" and its smooth compact perturbation 

Let us consider a stationary solution for a = 0, £ = x, choose a point xo > 0, a 
constant Ci = u(x + )p(x + ) > 0, and construct on the semi-axis x > x + a solution z + (x). 
Analogously for the semi- axis x < x_ we choose Ci = u(x)p(x^) < and construct a 
solution z-(x). Thus, outside of the segment [x_,x + ] we define a solution 

Z-(x), x G (—00, x_), 
z + [x), x G (x + , +00). 

For the sake of simplicity we set x_ = —x + , c\ = —k for z-{x) and c\ = k for 
z + (x), where k = const > 0. 

In their turn, the density, pressure and velocity can be found as 

p=*&., p = c 2 - k j, u=*^p. (3.5) 



Formation of singularities in solutions to ideal hydrodynamics of inelastic gases 



9 



It is very attractive to choose x + = x_ = and to construct a piecewise continuous 
solution like a solution of a "Riemmann problem" with non-constant left and right 
states. Nevertheless, it can be readily shown that the Hugoniot conditions do not hold 
on the jump. Indeed, the components p and p are continuous at the point x = 0, having 
a jump in derivative, however, the velocity itself has a jump [u] = -Mr > w^y (see 
Figs.l - 3). The value of [u] tends to zero as k — > oo, however the Hugoniot condition 
[pu] = does not implement in the origin x — 0. Of course, we can choose C\ such that 
the density and pressure have jumps in the origin and consider a^O, nevertheless the 
careful analysis shows that the Hugoniot conditions do not hold anyway. 




Figure 1. The density 




Figure 2. The pressure 



Therefore we choose R(0) > x + , and define smooth initial data such that for 
\x\ > -R(O) they coincide with p(x), u(x), p(x) and on the segment |x| < -ff(O) they 
are arbitrary smooth functions p(0, x) > 0, p(0, x) > 0, u(0, x). We will call this type of 
initial data the compact perturbation of nontrivial finally steady state. Let us note that 
for k = we get the trivial zero-state solution. 
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Figure 3. The velocity 
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5.5. Breakdown of the compact perturbation of the nontrivial "finally steady state" 

Let us denote the perturbed region by B{t) := {x\ \x\ < R(t)}, and consider the analogs 
of functionals used in the previous section: 



p(t,x)(u(t,x), x) dx, 



G(t) = - J p(t,x)\x\ 2 dx, F(t) 



B(t) 
1 



B(t) 



p(t,x)\u(t,x)\ 2 dx, M(t) 



B(t) 



E k {t) 



P(t)= / p(t,x)u(t,x)dx, S(t) = 

B(t) B{t) 

The following theorem holds: 



p(t,x) dx. 

B(t) 

K(t,x)p(t,x) dx. 



Theorem 2 Let P 2 (0) > 



SMc 2 
k ' 



n = 1. Then there exists no globally in t smooth 



perturbation of the nontrivial steady state for system ( fi. 1\) - $1.3\) . 



Proof. First of all we note that system (II. ip - (II. 3p is hyperbolic and therefore the 
speed of boundary of the perturbations is equals to |ii| +V S , where V s = y/pp, the sound 
speed. As follows from (13. 5p . 



T/ k <r ° 2 

V s = -< — := a, 

p K 

where we use the estimate p > — . 

r — C 2 

Then we have 

F'{t) = 2E k + [ (p(t, x) - p{t, R(t))) dx, 

M'{t) = 0, P'{t) = 0. 



(3.6) 



(3.7) 



B(t) 
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Further, the Holder inequality implies 

F 2 (t) 
E k (t) > =-M, (3.8) 

W_ 4G(t) V ^ 

and one can estimate 

G{t)< l -R\t)M{t)< l -R\t)M^). (3.9) 

Further, the Jensen inequality yields 

p(t, x) dx > (2R(t)) 1 -' / KX X (S(t)y. (3.10) 

B(t) 

As in the proof of Lemma [5] we can show that for sufficiently large t the function 
S(t) > S_ = const > 0. Then we notice that (13. 3p implies 

Ak 2 r 2 
z \ x ) ~ — To~ x } a; — )• oo. (3.11) 

Therefore taking into account (13. 5p we have 

/ p(t,R(t)))dx < % f z(R(t))dx ~ %lBT x (t), t->oo. (3.12) 
J k J A 2 

B{t) Bit) 

Thus, from (13. 6p - (13. 12[) and the estimate 

m) > p{t) - Pm 



2M(t) 2M(0) 

we have 

~, - F 2 (t) Kl~ l SI Ac 2 

{t) ~ k[t) + 4 M (0) (R(0) + at) 2 + (R(0) + crt)^ A 2 (R(0) + at) (3 ' 13) 

- 2M(0) 4M(0)( J R(0) + cxt) 2 
where A(t) — V 0, t — V oo. 

As we can see from (13.131) . beginning from a certain t > the function F(t) > 
4 M(0) -R(t ) cr, moreover, 

F 2 rti 

F'(t) > -^ — , t>t . (3.14) 

V ; ~ 4M{0){R{0) + at) 2 V 7 

Integrating (13. 14j) from t = to we get 

4M(0)F(t )i?(to)(i?(to) + <rt) 



F(t) > 



4M(0)i? 2 (t ) + (4M(0)i?(t )a - F(t ))t 

Thus, F(t) blow-ups at a finite time. This contradicts to the inequality F 2 (t) < 
4G(t)E(0). The theorem is proved. □ 

Remark 3 The idea of the method is due to [10], where it was proved that the compact 
smooth perturbation of a constant state of gas dynamics equation can not be globally 
smooth in time. 
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4. Exact solution with singularity 

Naturally there arises a question on a type of predicted singularity. In particular, in the 
remarkable papers [3J, [I] for the one- dimensional case the authors employ Lagrangian 
coordinates and derive a broad family of exact non- stationary non-self-similar solutions. 
These solutions exhibit a singularity, where the density blowups in a finite time when 
starting from smooth initial conditions. Moreover, the velocity gradient also blowups 
while the velocity itself and develop a cusp discontinuity (rather then a shock) at the 
point of singularity. This approach is partially extended to the 2D case in [5] . 

Here for any spatial dimensions we construct a simple family of solutions to the 
system (II. ip - (11. 3p having a singularity in the density whereas other components are 
continuous. Indeed, if we substitute in (II. ip - ( II. 3p 

u(t,x) = a(t)x, pit, x) = (5{t) \x\ q , p(t, x) = s(t) \x\ l , (4.1) 

where x is a radius- vector of point, we obtain 

q = -1, 1 = 1, P{t) = P = const > 0, (4.2) 

and ait), s{t) > satisfy the following system of nonlinear ODE: 

a'(t) + a 2 (t) + ^- = 0, (4.3) 

Po 

s'(t) + (7 + 1) n sit) ait) = -A (3 l/2 (t) s 3/2 (t). (4.4) 

This system has a unique equilibrium (a(i) = 0, sit) = 0), it is unstable. One of 
its solutions is very simple: s(i) = (p = 0), a(t) = (t + a~ 1 (0))~ 1 . An analysis of the 
phase portrait shows that if s(0) > 0, then a(i) — > — oo, s(i) — > +oo for all a(0). 

Let us prove this fact in a different way. We consider a symmetric material volume 
B(t) containing the origin x = and use the denotation of Sec l3.31 We can see that in 
spite of the singularity in the component of density all integrals below exist. Thus, due 
to the structure of solution (14.11) . (14.21) we have 

&{t) = F{t) = 2a{t)G{t), (4.5) 

F'(t) = 2 E k (t) + / ipit, x) - pit, R(t))) dx<2 E k (t) = 2 a 2 (t) G{t). (4.6) 

B{t) 

As follows from (14. 5p . (14. 6p . the velocity gradient obeys the inequality 

a\t) < -a 2 it), 

therefore a(t) < (t + a~ l (to))~ 1 , and in the case a (to) < we can see that a(t) — > — oo 
as t -y -a _1 (*o)- Since (@3D - (B3D result that a'(t) < -s(t)/f3 < -e < (the latter 
inequality follows from the uniqueness theorem), a(t) become negative in a finite time. 
The proof of over. 

We see that in the presence of a stationary singularity in the component of density, 
a balance between velocity and pressure arise. Generically the velocity and pressure 
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blow up in a finite time. Thus, the components of this "black hole" solution solution 
can collapse in different moments of time. The singularity of density in the origin is 
integrable for the dimension n > 2. 

Let us remark that for usual gas dynamics the solutions of such kind with a "linear 
profile" of velocity are well investigated (e.g. [9], Chapter IV, Sec. 15). 

The author is indebted to B.Meerson for attracting attention to the problem and 
thanks N.Leontiev and V.Shelkovich for a helpful discussion. 
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